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The behavior of liquid water under an electric field is a crucial phenomenon in science and engi- 
neering. However, its detailed description at a microscopic level is difficult to achieve experimentally. 
Here we report on the first ab initio molecular-dynamics study on water under an electric field. We 
observe that the hydrogen-bond length and the molecular orientation are significantly modified at 
low-to-moderate field intensities. Fields beyond a threshold of about 0.35 V/A are able to dissociate 
molecules and sustain an ionic current via a series of correlated proton jumps. Upon applying even 
more intense fields (~ 1.0 V/A), a 15 — 20% fraction of molecules are instantaneously dissociated 
and the resulting ionic flow yields a conductance of about 7.8 • cm -1 , in good agreement with 
experimental values. This result paves the way to quantum-accurate microscopic studies of the 
effect of electric fields on aqueous solutions and, thus, to massive applications of ab initio molecular 
dynamics in neurobiology, electrochemistry and hydrogen economy. 



Water is the most abundant compound on Earth and 
plays a fundamental role in the realm of natural and 
life sciences, at the interface between a variety of dis- 
ciplines such as physics, chemistry and biology [H [5], 
with a seemingly infinite range of technological applica- 
tions. Despite the interest and the enormous amount 
of fundamental and applied research work, a complete 
understanding of its unique behavior is still lacking [5J. 
The peculiar properties of water are due to its molecular 
structure and dipole moment and, as a consequence, to 
its hydrogen-bonding properties [7] . Among the physical- 
chemical phenomena of greater interest, particularly rel- 
evant is the phenomenon of autoprotolysis, i.e. the dis- 
sociation of a water molecule according to the following 
reaction [5]: 

2H 2 # OH~ + H 3 + (1) 

where two molecules form a pair of hydronium and hy- 
droxide ions which then separate, thus determining the 
pH of the substance. This effect is of paramount impor- 
tance because it is a key step in all processes in which 
water is involved as more than just a passive solvent. The 
kinetics of the autoionization process has been studied in 
some experiments [9] , but it is an extremely rare reaction 
which takes place on the femtosecond time scale of molec- 
ular motions |10j . thus making its microscopic investi- 
gation very difficult. On the other hand, a deeper un- 
derstanding of this reaction is crucially important in the 
presence of an external electric field, with consequences in 
disparate domains, from neurobiology [TT] to electrolytic 
batteries [H] and hydrogen-based technology [15] . For 
example, the activity of neurons is governed by the ex- 
change of ions through the cell membrane, which is driven 
by the potential difference between the interior of the 
cell and its outside environment. The first experiments 



of forward and backward conversion of chemical energy 
into electric energy date back to 1789; the first battery 
was invented by Volta in 1800 and then used to carry out, 
for the first time, the electrolysis of water. Since then, 
electrochemistry has become a stand-alone branch of sci- 
ence and technology, with a huge impact on everyday life 
and more to come on the nascent hydrogen economy. 

The chemical processes governing the autoprotolysis of 
water in the absence of an external field are in principle 
within the reach of ab initio molecular-dynamics simula- 
tions. However, because of the large timescales involved, 
a straightforward molecular-dynamics approach turns 
out to be rather ineffective. The problem was overcome 
using a combination of Car-Parrinello (CP) molecular dy- 
namics |14j and transition-path sampling |15j . which con- 
firmed that water dissociation is an extremely unfavored 
event from the energetic point of view. In a milestone 
study Geissler et al. showed that the rare autoprotolysis 
of a water molecule is induced by occasional large fluctu- 
ations of the electric field generated by the surrounding 
solvent, i.e. by other water molecules [TBJ. To test their 
theory, these authors adopted a clever scheme in which 
the occurrence of such rare fluctuations was favored by 
ad hoc semi-random moves of protons along hydrogen 
bonds; they showed that this method is equivalent to 
applying an external electric field and proved that such 
a field can dramatically change the energetic balance of 
the dissociation reaction, favoring the occurrence of the 
latter. The nascent hydronium and hydroxide ions then 
separate along a chain of hydrogen bonds via Grotthuss 
transfer events [T7]. It also turns out that their recom- 
bination is made possible through a collective compres- 
sion of the water wire bridging the ions |18) . However, 
when the hydrogen-bond chain still connecting the two 
ions breaks, their rapid recombination is inhibited and 
from such a state ions can move further apart and sep- 



2 



arate completely. Based on this result, recent classical 
molecular-dynamics simulation studies have shown that 
the important contributions to the electric field acting on 
OH bonds stem from water molecules less than 7 A away, 
corresponding to charges on water molecules of the first 
and second hydration shell [TS], and that they are able 
to generate instantaneous fields up to 2.5 — 3.0 V/A. 

Despite the major progress that such works represent 
in this domain, it has not been possible to establish, so 
far, a quantitative relation between the degree of water 
dissociation and the intensity of the electric field. Sev- 
eral studies based on the classical molecular dynamics of 
water molecules under an external electric field are al- 
ready present in the literature, describing, for example, 
the effect of an electric field on the phase diagram [20] 
or the phenomenon of electrofreezing [3TJ [25]. How- 
ever, such studies are typically based on rigid-molecule 
non-polarizable potential models, which are unable to 
describe the dissociation of water molecules. On the 
other hand, the effect of an electric field in ab initio 
Density Functional Theory (DFT) codes is traditionally 
taken into account through a sawtooth potential that is 
quite satisfactory in low-dimensional systems but which 
cannot be used in bulk systems. We present here the 
first, to our knowledge, ab initio study of the properties 
of bulk water in the presence of intense external elec- 
tric fields. The study was carried out using Quantum 
Espresso |23j . a well-known open-source package for the 
DFT calculation of electronic-structure properties, which 
uses a plane- wave basis set and pseudopotentials. Macro- 
scopic polarization and finite electric fields are treated by 
this code with the modern theory of polarization (Berry 
phases) [23] . Previous ab initio calculations of water in 
an external field were typically carried out on small clus- 
ters in the vacuum [25 . All the technical details of the 
calculation are reported in the Supplementary Material 
file [2B]. 

In order to test the reliability of the method, we first 
checked the effect of a uniform electric field by calcu- 
lating the field-induced electronic polarization at fixed 
atomic positions for several water configurations as in 
Ref. [24]; the estimate we obtained (£ m s» 1.74) is in 
good agreement with the experimental data as well as 
with previous theoretical results (we refer the reader to 
Ref. [25] for more detailed information). We also calcu- 
lated the structural properties of the system at zero field, 
such as the radial distribution functions, and found that 
they reproduce with fair accuracy those reported in the 
literature [2B]. As shown in Ref. [37], the evolution of 
the structure of water and of its H-bonding properties 
are better described by the instantaneous distribution 
of the distance between an oxygen atom and its third 
hydrogen neighbor, i.e. typically the closest H neigh- 
bor covalently bonded to a nearby oxygen atom. Wc 
report such a distribution, for values of the external elec- 
tric field up to 1 V/A, in the top panel of Fig. [I] to- 
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FIG. 1: Distributions of istantaneous O — O and O — H dis- 
tances for electric fields in the range between zero (purple 
curves) and 1 V/A (red curves). The insets show the evolution 
of most probable distances with the field intensity. Top panel: 
intermolecular distance between an oxygen atom and its third 
hydrogen neighbor; central panel: intramolecular distance be- 
tween an oxygen atom and its second hydrogen neighbor; bot- 
tom panel: oxygen-oxygen closest-neighbor distance. The ex- 
perimental data were smoothed with the Savitzky-Golay filter 
method, as typically shown for the zero-field curve in the bot- 
tom panel. 



gether with the intramolecular O — H distance distribu- 
tion of the second hydrogen neighbor (central panel) , and 
the — closest-neighbor distance distribution (bot- 
tom panel). The effects produced by increasing elec- 
tric fields on the first distribution (top panel) are a sys- 
tematic broadening and a significant shortening of the 
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FIG. 2: Counts of OH and H30 + ionic species with lifetimes 
as abscissa, for field intensities 0.25 V/A (yellow), 0.35 V/A 
(blue), and 1 V/A (red). 

H-bond for fields up to about 0.35 V/A; for larger in- 
tensities, the maximum of the distribution weakly shifts 
again to larger and larger values while a secondary peak 
emerges at about 1.1 A, which can be associated to a co- 
valently bonded hydrogen in a H30 + ionic species. The 
effect of the field is also manifest on the intramolecu- 
lar covalent bond between an oxygen atom and its sec- 
ond hydrogen neighbor, as shown in the central panel 
of Fig. [T] the associated distribution becomes more and 
more delocalized while the most probable distance ex- 
hibits a nonmonotonic trend, opposite to that observed 
above for the hydrogen bond. Note also, in the bottom 
panel of Fig. [I] the monotonic shortening of the O — O 
distance up to an amount of about 0.14 A, a shift com- 
parable to that observed in water under a pressure of 
about 20 GPa 28J. In classical simulations intermolec- 
ular distances are unaffected while an increased number 
of H-bonds has been observed with longer lifetimes [21] . 
We interpret this partial discrepancy as a clear indication 
of a larger, field-induced O — H intramolecular polariza- 
tion, which strengthens the intermolecular O — H inter- 
action and which cannot be accounted for in the classical 
models most commonly used to describe water. We also 
observed, indeed rather unsurprisingly as this effect has 
been already predicted by classical simulations [3T] , that 
water molecules tend to align their dipoles along a direc- 
tion parallel to the field [2"6] . 



FIG. 3: Snapshot of a typical dissociation-diffusion mech- 
anism; the molecules involved in the process are rendered 
with red coloring, while the surrounding molecules have gray- 
colored oxygen atoms. In both cases, hydrogens are de- 
picted in white. In particular, a HaO + — H2O — OH - entity 
is clearly visible and its relevant hydrogen-bonds are indicated 
with a dashed white line. The electric field is oriented along 
the positive z direction. 



According to the present calculations, water molecules 
start to dissociate for field intensities around 0.25 V/A. 
Protons participating in H-bonds, mostly aligned with 
the external field, start to jump occasionally back 
and forth along the bond, thus creating instantaneous 
H 3 + — OH~ pairs. Such events are still rare and pro- 
duce extremely short-lived (< 10/s) ionic pairs, the field 
being not strong enough to permanently sustain instan- 
taneous ionic dipoles for times long enough to start a real 
proton diffusion (see Fig. [2]). Only for fields of about 0.35 
v/A we started to observe much more frequent molecular 
dissociations, suggesting that the balance of the reaction 
described in Eq. [I] more and more favors the products in 
the right term at the expense of the reactants in the left 
term, i.e. the neutral state. We note that the above field 
threshold is in very good agreement with experiments. 
In a work of a few years ago [3Q|, water dissociation was 
observed to occur under external fields of 0.32 and 0.44 
V/A, the threshold depending on the temperature. A 
very recent experimental study provided a thorough re- 
view of the process of water ionization in high electric 
fields induced by Pt electrodes [3T]. In particular, an in- 
teresting analysis of the interplay between the local field, 
as estimated by Geissler et al. [TB] , and the external field 
concludes that the latter should enhance water ionization 
for intensities of 0.3 — 0.6 V/A, in excellent agreement 
with our results. 

Once the external field is strong enough to induce a 
significant dissociation, ionic lifetimes become dramati- 
cally longer, ranging from a few femtoseconds to more 
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than 0.1 ps for both H30 + and OH~, as reported in 
Fig. [2] Indeed, at this point the dissociation reactions 
are typically followed by a series of proton jumps from 
the newly formed hydronium to the adjacent H-bonded 
neutral molecule along the field direction, forming a new 
hydronium from which a (different) proton can jump, and 
so on and so forth. In other words, ionic diffusion (and 
thus conduction) is activated by the dissociation reac- 
tion described in Eq. 0, but the presence of the field 
breaks the symmetry of the physical space, providing a 
spatial direction for the relevant proton hopping events, 
and thus favoring (or disfavoring) their respective reac- 
tion equilibria according to the relative physical position 
of the reactants, with respect to each other and to the 
field. This effect can be visually pictured as follows: 



H 2 -I- H 2 — > OH" + H 3 0^ 



(2) 



where the right-directed oriented field induces a proton 
jump from the left water molecule to the right one. By 
analogy beyond the dissociation threshold the recombi- 
nation processes more likely occurs as: 



H 3 0+ + OH" 



H 2 + H 2 0, 



(3) 



i.e. the extra proton comes from an hydronium in the 
opposite direction of the hydroxide with respect to the 
field. Besides these two steps, ionic conduction is sus- 
tained by other processes which can be decomposed in 
terms of the following elementary proton jumps: 



H 3 0+ + H 2 — > H 2 + H 3 H 



and 



H 2 + OH~ 



OH" 



H 2 0, 



(4) 



(5) 



where, again, proton hopping along the field direction 
is privileged. We note that proton jumps are signifi- 
cantly correlated; in particular, protons originally form- 
ing donor H-bonds and belonging to adjacent molecules 
jump almost simultaneously as their respective acceptor 
molecule, in a Zundel-to-Zundel mechanism very similar 
to the one described in ab-initio simulations [16] as well as 
in experiments [29] . The mechanism of correlated proton 
jump diffusion is illustrated by a snapshot of our simula- 
tion cell in Fig. [3j Dissociated molecules recombine with 
protons diffusing from neighboring molecules, thus con- 
tributing to the global liquid conductivity. We observed 
that even for fields of 1.0 V/A adjacent H 3 0+ - OH" 
pairs are not viable, unless they are separated by an H 2 
molecule, and that the Zundel-to-Zundel diffusion mech- 
anism is dominant. Such correlated processes can involve 
even more than two protons at the same time, but still 
the decomposition of the diffusion mechanism in the in- 
dividual elementary steps of Eqs.pflBl is possible, and 



allows a more quantitative description of the whole phe- 
nomenon and of the corresponding rates. To this end, 
we reported in Table [T] the number of such elementary 
events per molecule and per picosecond. We note that at 
the dissociation threshold (0.35 V/A) a single molecule 
undergoes, on average, about 3 proton-hopping events 
in a timespan of 1 picosecond, about 6 at 0.5 V/A, up 
to about 10 at the highest field considered in this work, 
and that the fraction of instantaneously ionized molecules 
AiV increases accordingly. In particular, at the most in- 
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0.148 
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0.067 


0.108 


0.0046 
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0.952 


0.872 


0.743 
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0.0408 
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1.049 


1.808 


1.856 


0.0828 


1.00 


2.341 


2.082 


2.502 


2.728 


0.1556 



TABLE I: Calculated event rates and ion fractions for dif- 
ferent field intensities, per molecule and per picosecond of 
dynamical evolution: dissociation (Eq.pl); recombination 
(Eq.Q); hydronium-to- water hopping (Eq.Q); water-to- 
hydroxide hopping (Eq. ([5J ) . The average fraction of ionized 
molecules is reported in the last column. 



tense field a fraction AN ss 0.15 — 0.20 of the total num- 
ber of molecules are instantaneously ionic, as appears 
from the rightmost column of Table [I] giving rise to a 
ohmic conductivity of about 7.8 f2 _1 • cm" 1 ; this lat- 
ter finding is in fairly good agreement with the available 
experimental values [32H34] of water dissociation under- 
pressures larger than 20 GPa, as analogously observed 
in ab initio CP calculations carried out in similar ex- 
treme conditions [551 123 . In particular, the authors of 
Refs. 375] and [21! observed in shock-wave experiments 
a conductivity of about 30 ft -1 ■ cm -1 , which they at- 
tributed to a massive presence of ionic species; they also 
estimated that fully ionized water (i.e. 100 % of disso- 
ciated molecules) should have, at those regimes, a con- 
ductivity of about 36 O -1 • cm -1 . Our result thus seems 
to indicate that strong electric fields, yet similar to those 
instantaneously generated in neutral water by the sur- 
rounding molecules [19] , are able to generate ionic species 
and to induce an ionic current whose intensity is compa- 
rable to that observed in different experimental condi- 
tions and setups. The calculated ionic current- volt age 
characteristic of water is finally presented in Fig. |4j 

We conclude that the ab initio Car-Parrinello approach 
is able to catch the essential features of molecular disso- 
ciation and diffusion mechanisms of bulk water under the 
action of an intense electric field. The results illustrated 
in this Letter show that the method employed not only ef- 
ficaciously describes the dissociation of the molecule but 
further accounts for the presence of the forming ionic 
species whose effect, even if short lived, is very impor- 
tant in the dynamics of ions in solution. It appears that 
ab initio simulations of aqueous/molecular systems under 
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FIG. 4: Ionic current-voltage characteristic for a cubic cell of 
12.42 A side. The dots represent the calculated values, while 
the solid line is a guide for the eye. 



an electric field are now feasible with a numerical accu- 
racy which may already allow for a comparison with the 
available experimental data. A stimulating scenario of 
new investigations opens up with a range of applications 
of potentially enormous reach, from planetary science, 
to biology, electrochemistry and the fast-growing field of 
hydrogen-based economy. 
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